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For use in the computation of thermodynamic functions, analytical descriptions are given for the 
following properties: the orthobaric densities and saturation temperatures; the heats of vaporization; 
the specific heats of saturated liquid; and the thermodynamic functions for ideal gas states. 
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List of Symbols 

Subscripts c and t refer to critical and liquid triple 
points. 

Subscripts g and / refer to vapor and liquid in 
equilibrium. 

a, b, c, constant coefficients 

A; constant coefficients 

/3, an exponent, 0.36 

C r , molal heat capacity at constant volume 

C ; >, molal heat capacity at constant pressure 

C (r , molal heat capacity of saturated liquid 

d, molal density 

A, any relative deviation, (y/y ca ic-0 

E, the internal energy 

e, exponent in various expressions 
//, the enthalpy 

J, the joule 

A:, conversion factor, 101.325 .///-atm) 

/, the liter, 1000 cm 3 

mol, 16.042 grams of methane 

P, pressure, 1 atm = 101325N/m 2 

3>(r) , (E° -£§)/RT, internal energy function 

Q r , heat absorbed in a reversible change of state 

R, the gas constant, 8.3143 J/(mol-K), [26] 

p, d/df, density reduced at dt 

S, the entropy 

cr, djdc, density reduced at d c 

T, temperature, K, IPTS (1968) 

r. s (p), liquid- vapor coexistence temperature 

r, T/Tc, temperature reduced at T ( 

u, (T c — T)l(T r — Tt), normalized temperature, 
origin at T c 



v, \jd, molal volume 

w, (T c IT—l)l(TclTt—l), inverse normal temp., 

origin at T ( 

oj, e/x, an argument in eq (8) 

x, 77400, reduced temperature for eq (8) 

F, function of C (r , defined on figure 5 

z, (1 — r), reduced temperature with origin at T ( . 

The number of digits, given for coefficients in the 
empirical functions of this report, is chosen to yield 
calculated results having at least the same number of 
figures as exist in the data to be described. 

1 . Introduction 

The current need for thermodynamic properties of 
methane was described in our preceding report [I]. 1 
The present report concludes our selection of pub- 
lished properties data for coexisting vapor and liquid, 
and includes also a description of the ideal gas thermo- 
functions. All of these properties are placed in analyti- 
cal form directly useful for thermal computations. 

Temperatures are on the IPTS (1968), as in [1], and 
we use the same fixed-point constants: 7^ = 90.66 K, 
P, = 0.1151 atm, (0.01166MN/m 2 ), d*= 28.0536 mol// 
(liquid); and T c = 190.53 K, P c = 45.346 atm, (4.5947 
MN/m 2 ), r/ r = 10.15 mol//. 

2. The Orthobaric Densities 

For the densities of coexisting vapor and liquid 
methane the preponderance of published data is on 
the liquid phase [2, 3, 4, 5, 6, 7, 8]. For the vapor 
phase, therefore, we have estimated data at densities 
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1 Figures in brackets indicate the literature references at the end of this paper. 
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up to 1/3 of dc by use of the vapor pressure equation 
and the virial equation of state. We used the repre- 
sentation of second virial coefficients based on data of 
Sengers, Klein, and Gallagher, as well as the third 
virial coefficient, both as described in [1]. 

Analytical description of the vapor densities is given 
by our new form, constrained to qualitatively accepta- 
ble behavior near the critical point [9]. To facilitate 
selection of significant figures in the coefficients, we 
use variables u and w normalized to a maximum value 
of unity (at the triple point: see list of symbols). Vapor 
densities now are described by use of exponent 
0=0.36, 



In (p g lpc)=A 1 -w + A 2 -uP + ^Ai-u i - 

i=3 



(1) 



A l =- 


-27.42 


792 


A 5 = - 


-10.44 


279 


A 2 = 


-1.56 


983 


A 6 = 


23.44 


911 


A 3 = 


10.16 


428 


A 7 = - 


-18.38 


419 


A 4 = 


10.16 


575 


A s = 


7.56 


705 



The value for exponent ft was selected by consider- 
ing the published values, (3 = 0.367 [10] , and £ = 0.3566 
[11]. The number of terms in (1) was selected to bring 
deviations down to the level of estimated uncertainty 
in data derived from the vapor pressure and virial 
equations. Individual relative deviations of vapor 
densities from eq (1) are plotted on figure 1. The rms 
deviation for 55 points is 0.21 percent. 

Saturated liquid densities are described by our 
simple form [9] using argument z = ( 1 — t) , 

pilpc-l = a-z+6-^+c- exp (-e-r 2 /z), (2) 
a = 0.62 253 c = -0.08 367 
6 = 1.89 542 e= 0.6 

Inconsistent data were eliminated successively, leav- 
ing 60 points with an rms relative deviation of 0.24 
percent. Figure 2 gives individual relative deviations 



from eq (2), and identifies the sources. At the triple 
point 90.66 K, eq (2) gives a liquid density d t = 28.0536 
mol// resting primarily on the data of [8] , as seen on 
the figure. 

Figure 3 gives the coexisting densities of vapor and 
liquid methane along the abcissa, calculated by eqs 
(1, 2). The ordinate gives the saturation temperature, 
T s (p). This plot provides a quick, rough estimate of 
the data. Conspicuously absent is the line for densi- 
ties of freezing liquid. Some values recently have been 
estimated to 94.5 K (150 atm), including the volume 
change and heats of melting [12]. 

The rectilinear diameter may be derived by use of 
eqs (1, 2) for interpolation of d g and di to the same tem- 
peratures. For 39 points at 112 ^ T^ 188 K, our rep- 
resentation in mol//, 

(dg + di)/2 = a + b-z+c exp (- e/z) , (3) 

a=10.156, 6 = 7.209, c = 0.183, 6 = 3/4, 

yields an rms deviation of 0.05 percent and a critical 
density d c = 10.156 + 0.02 mol//, as compared with 
d c — 10.15 adopted for present work. 

The relation, T s (p), is important in the nonanalytic 
equation of state [13] . For the combined data of figures 
1 and 2 we obtain coefficients for the expression, 



7yr s =i+|o— i| 3 -F( P ), 

F(p)^A 1 -ln(llp)+f i Arp i -\ 



(4) 



A x = 


0.19 3073 


A e = - 


-112.50 780 


A 2 = 


-0.34 282 


A 7 = 


141.02 284 


A 3 = 


2.44 483 


A s = 


-89.60 804 


A 4 = - 


- 12.24 696 


A 9 = 


22.74 904 


A,= 


48.68 975 


rms = 


0.20 % 



The signs and magnitudes of these coefficients are 
comparable with values found for other substances [9], 
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FIGURE 1. Relative deviations of the densities of saturated methane vapor from eq (1). 

Open circles from the virial surface and vapor pressures; filled circles from [4]; circles with tails from [5]. 
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FIGURE 2. Relative deviations of the densities of saturated liquid methane from eq (2). 

Filled circles from [8]; open circles from [7]; tails up from [5|; tails down from [4]. 
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FIGURE 3. Methane saturation temperatures as a function of vapor and liquid densities up to the 

triple point. 



3. The Heats of Vaporization 

Some data for the heats of vaporization, A// r , have 
been reported [14, 15]. We now use the Clapeyron 
equation to calculate these heats by use of our vapor 
pressure equation [1] and eqs (1 and 2) of the present 
report, 

&H v =k-T-(dPldT)-(vg-vi). (5) 

Column 1 of table 1 gives the same rounded tempera- 
tures used in [14, 15]. (Experimental uncertainties do 
not demand conversion to the IPTS-1968.) Second and 
third columns of table 1 give the published data in 
kilo-Joules/mole. The fourth column gives our results 
via (5). The remaining three columns are the derived 
data used on the right of (5). Figure 4 gives the outline 
for our calculated results. 

For some computations it would be convenient to 
use a more direct description of AH r than given by (5). 
For the 21 points in the fourth column of table 1 the fol- 
lowing representation gives a maximum relative devia- 
tion A = — 0.04 percent at 120 K; deviations of 0.00 



percent at T= 185 and T= 190 K; and an rms devia- 
tion of 0.02 percent, 

bHv^JZAcu*!*, kj/mol, (6) 

i=l 



At = 


7.88 860 


A 4 =- 


- 152.63 520 


A 2 = - 


-12.94 022 


A b = 


136.57 252 


A s = 


76.33 098 


A K = 


-46.55 421 



4. Specific Heats of Saturated Liquid 

Specific heats of the liquid phase, along the co- 
existence path, are useful for computing thermo- 
dynamic properties into compressed liquid states. 
Data for Co- are derived from observations on the two 
phase, liquid-vapor system at constant volume by use 
of accurate PVT data for the two phases [17]. Pub- 
lished data for methane are quite uncertain because 
the authors failed to give either their experimental ob- 
servations or a quantitative description of their der- 
ivations [14, 15]. The data are shown by figure 5 in 
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Table 1. Heats of vaporization for methane, kj/mol 



7\ K 


A//, 


dPfdT 

atm/K 


Vg 
//mol 


Vl 


[14] 


[15] 


Calc 


//mol 


90.66 






8.661 


0.0147 


64.1399 


0.0356 


95.00 






8.597 


.0227 


39.4535 


.0361 


100.00 


8.571 


8.494 


8.506 


.0353 


23.8300 


.0367 


105.00 


8.433 


8.355 


8.394 


.0522 


15.1377 


.0373 


110.00 


8.300 


8.213 


8.260 


.0741 


10.0373 


.0379 


115.00 


8.162 


8.063 


8.110 


.1014 


6.9043 


.0385 


120.00 


8.025 


7.899 


7.943 


.1344 


4.9011 


.0392 


125.00 


7.878 


7.724 


7.761 


.1734 


3.5743 


.0400 


130.00 


7.731 


7.531 


7.563 


.2186 


2.6672 


.0408 


135.00 


7.565 


7.326 


7.347 


.2702 


2.0295 


.0416 


140.00 


7.394 


7.104 


7.112 


.3283 


1.5698 


.0425 


145.00 


7.203 


6.862 


6.856 


.3930 


1.2310 


.0436 


150.00 


6.987 


6.598 


6.578 


.4645 


0.9763 


.0447 


155.00 


6.740 


6.305 


6.273 


.5431 


.7814 


.0460 


160.00 


6.468 


5.987 


5.937 


.6292 


.6296 


.0475 


165.00 


6.153 


5.619 


5.563 


.7232 


.5093 


.0493 


170.00 


5.784 


5.197 


5.134 


.8262 


.4121 


.0514 


175.00 


5.339 


4.699 


4.623 


.9397 


.3316 


.0542 


180.00 


4.776 


4.021 


3.980 


1.0664 


.2625 


.0579 


185.00 


3.933 


3.109 


3.088 


1.2122 


.1995 


.0636 


190.00 




0.782 


1.267 


1.4059 


.1269 


.0801 




FIGURE 4. Outline for the heats of vaporization of methane, kilo-Joules/ mol, derived via the Clapeyron 
equation with analytical descriptions of the physical properties. 



coordinates C (T 'U € versus T, the weighted form used for 
least squares determination of coefficients in 

C (T = a + b'U + clu% J/(mol K), (7) 

a = 29.250 c= 18.924 
b= 5.560 6= 0.55 

The value for exponent e was selected by analogy: our 
oxygen data yield e = 0.54 by trial, and our fluorine data 
give 6 = 0.59 [18]. Present methane data yield 
e=1.0±0.1. The line on figure 5 is calculated from 



(7). The rms relative deviation is 2.6 percent for 33 
values of C lT . 



5. The Ideal Gas States 

Thermodynamic functions for hypothetical ideal gas 
states at one atmosphere pressure may be derived for 
simple molecules as a function of temperature from 
spectroscopic data by statistical methods [19, 20]. For 
our thermodynamic computations with PVT data, 
these ideal gas functions are invaluable because their 
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Figure 5. The function Y = C, r *u 0M [J I moldeg) for specific heat of liquid methane along the saturation 

path, using u = (T c — T)/(T C — T t ). 



Filled circles arc from [15]; op< 



cles Iron. [14]; circles with tails from [16]. The 1 



accuracy exceeds that which can be attained by experi- 
mental specific heat measurements. 

Data for methane first were derived in API Project 
44 [21], and work to 1961 has been reviewed [22]. For 
present work we select the recent results of McDowell 
and Kruse [23], and give them an analytical repre- 
sentation in the range 60^7^400 K. With conven- 
tional notation [19], abbreviate the internal energy 
function by *(f) « (E -E$)/RT, and define the 
argument x — 7/400 with maximum value of about 
unity for our range of T. 

A formulation with six constants was used for nitro- 
gen [24]. For methane [23] and for oxygen [25] it does 
not give a highly precise representation of the specific 
heats. A more precise description is obtained with the 
power series, 

On the other hand, a modification of the form given in 
[24] diminishes specific heat rms deviations by a 
factor of about 1/2. Maximum deviations in C%/R 
(given below) become 0.02 percent in our present 
range of interest for methane, 80 ^T^ 360 K. This 
accuracy is at least an order of magnitude better than 
can be obtained from PVT data when using the thermo- 
dynamic computation, 



C v (p, T) = C°(T)-T- (d*PldT*)-dplp* 



The modified form (8) uses an argument, cj = e/x 
Constant e is found by trial, 



®(T) =Ax + Az-x^+Aa-x*!* +A 4 'X + As-a>l (e»-l) i 



[16]. The line is calculated 


,vith eq (7). 




e = 5.022 880 


A 3 = - 


-1.847 272 


Ai = 2.599 898 


A 4 = 


0.821 122 


A>= 1.444 942 


A s = 


4.720 791 



The specific heat is obtained via C{J = dE°/dT, 

C%IR =A X + l-ArxM+i • A* -x 2 ' 3 + 2-A 4 -x 

+ ^5-[oi/(c«- 1) | 2 -e*\ 

The entropy is obtained via its definition, 



(9) 



AS - Q r /T= [C v -dTIT+ P-dv/T] . 

Jy, 

Introducing the published value of S°/R at T\ = 60 K 
yields, 

S°IR=A + ]n (7760)+ A -In (*) +4v4 2 -* 1/3 

+ $-A 3 'X*l* + 2-A4-x + A*'[<al(e»--l)-ln(l — e-»)] 

(10) 
with constant A = 18.852 484. 

Specific heat at constant pressure, and the enthalpy 
function now are simply 



C°/jR = 1 + C°//J, 

(H Q -E§)/RT=1 + ®(T). 



(11) 
(12) 



(8) 



Table 2 gives temperatures in the first column. The 
next three pairs of columns give data and calculated 
values for <&(r), for Cj>/#, and for S°/R. 

Independent computations of above data by Harri- 
son et al. [27] give almost identical results at 90 ^ T 
< 800 K, further increasing confidence in their 
accuracy. 
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Table 2. 


Ideal gas 


functions for methane at 1 atm 




T, K 


*(T) 


Calc 


cm 


Calc 


S«IR 


Calc 


60.00 


2.9693 


2.9693 


3.0014 


3.0007 


15.9337 


15.9337 


80.00 


2.9774 


2.9774 


3.0019 


3.0021 


17.0849 


17.0849 


90.67 


2.9803 


2.9803 


3.0021 


3.0023 


17.5859 


17.5859 


100.00 


2.9823 


2.9823 


3.0024 


3.0023 


17.9779 


17.9779 


111.67 


2.9844 


2.9844 


3.0027 


3.0024 


18.4197 


18.4197 


120.00 


2.9857 


2.9857 


3.0029 


3.0026 


18.7077 


18.7077 


140.00 


2.9882 


2.9882 


3.0042 


3.0039 


19.3248 


19.3248 


160.00 


2.9904 


2.9903 


3.0075 


3.0075 


19.8597 


19.8596 


180.00 


2.9927 


2.9927 


3.0153 


3.0157 


20.3321 


20.3320 


200.00 


2.9956 


2.9957 


3.0309 


3.0313 


20.7559 


20.7559 


220.00 


3.0000 


3.0000 


3.0571 


3.0574 


21.1412 


21.1412 


240.00 


3.0063 


3.0063 


3.0966 


3.0966 


21.4958 


21.4958 


260.00 


3.0152 


3.0152 


3.1506 


3.1503 


21.8257 


21.8258 


273.15 


3.0227 


3.0228 


3.1943 


3.1939 


22.0316 


22.0316 


280.00 


3.0272 


3.0272 


3.2196 


3.2191 


22.1357 


22.1358 


298.15 


3.0412 


3.0412 


3.2947 


3.2942 


22.4030 


22.4030 


300.00 


3.0428 


3.0427 


3.3030 


3.3025 


22.4296 


22.4296 


320.00 


3.0620 


3.0619 


3.3996 


•3.3994 


22.7104 


22.7103 


340.00 


3.0850 


3.0849 


3.5077 


3.5080 


22.9803 


22.9802 


360.00 


3.1117 


3.1117 


3.6257 


3.6264 


23.2412 


23.2412 


380.00 


3.1420 


3.1421 


3.7516 


3.7527 


23.4947 


23.4947 


400.00 


3.1758 


3.1759 


3.8839 


3.8847 


23.7417 


23.7418 


450.00 


3.2737 


3.2737 


4.2328 


4.2291 


24.3369 


24.3369 



6. Discussion 

This report gives analytical descriptions of some 
methane properties, using existing and estimated 
data, as a base for preliminary computations of thermo- 
dynamic functions over the wide range from gaseous 
to compressed liquid states at low temperatures. It 
reveals deficiencies in the quantity of data available 
for densities of saturated vapor; in the precision of 
data available for saturated liquid densities; and hence 
uncertainty, difficult to assess, in the derived heats of 
vaporization. Data available for specific heats of liquid 
at coexistence cannot be assessed, due to a lack of 
published information. For thermofunctions in ideal 
gas states, the high quality of results derived by 
McDowell and Kruse suggests that they will stand for 
a long time. 

This laboratory is undertaking an experimental 
program to obtain reliable and self-consistent measure- 
ments over a wide range of conditions and to compute 
thermodynamic properties based on their formulations. 
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